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WITH  HEAT  TRANSFER  AND  PRESSURE  GRADIENT 
By  Clarence  B.  Cohen  and  Eli  Reshotko 


SUMMARY 

Stewartson 's  transformation  is  applied  to  the  laminar  compressible 
houndary-layer  equations  and  the  requirement  of  similarity  is  introduced, 
resulting  in  a  set  of  ordinary  nonlinear  differential  equations  pre¬ 
viously  quoted  hy  Stewartson,  hut  unsolved.  The  requirements  of  the 
system  are:  Prandtl  number  of  1.0,  linear  viscosity-temperature  rela¬ 
tion  across  the  boundary  layer,  an  isothermal  surface,  and  the  particular 
distributions  of  free-stream  velocity  consistent  with  similar  solutions. 
This  system  admits  axial  pressure  gradients  of  arbitrary  magnitude,  heat 
flux  normal  to  the  surface,  and  arbitrary  Mach  numbers. 

The  system  of  differential  equations  is  transformed  to  an  integral 
system,  with  the  velocity  ratio  as  the  independent  variable.  For  this 
system,  solutions  are  found  for  pressure  gradients  varying  from  that 
causing  separation  to  the  infinitely  favorable  gradient  and  for  wall 
temperatures  from  absolute  zero  to  twice  the  free-stream  stagnation  tem¬ 
perature.  Some  solutions  for  separated  flows  are  also  presented. 

For  favorable  pressure  gradients,  the  solutions  are  unique.  For 
adverse  pressure  gradients,  where  the  solutions  are  not  unique,  two 
solutions  of  the  infinite  family  of  possible  solutions  are  identified  as 
essentially  viscid  at  the  outer  edge  of  the  boundary  layer  and  the  re¬ 
mainder  essentially  inviscid.  For  the  case  of  favorable  pressure  gradients 
with  heated  walls,  the  velocity  within  a  portion  of  the  boundary  layer  is 
shown  to  exceed  the  local  external  velocity.  The  variation  of  a  Reynolds 
analogy  parameter,  which  indicates  the  ratio  of  skin  friction  to  heat 
transfer,  is  from  zero  to  7.4  for  a  surface  of  temperature  twice  the 
free-stream  stagnation  temperature,  and  from  zero  to  2.8  for  a  surface 
held  at  absolute  zero  where  the  value  2  applies  to  a  flat  plate. 
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INTEOIdJCTION1 

Factors  that  affect  the  development  of  laminar  boundary  layers  are 
pressure  gradient,  Mach  number,  and  heat  transfer,  plus  the  properties 
of  the  fluid  under  consideration.  Since  mathematical  complexities  pre¬ 
clude  solutions  of  this  problem  in  a  completely  general  fashion,  the 
literature  consists  largely  of  solutions  treating  particular  combinations 
of  these  factors.  For  the  flow  of  an  ideal  gas  over  a  surface  without 
pressure  gradient,  the  remaining  factors  have  been  taken  into  account 
very  completely  by  Crocco  (ref.  2)  and  Chapman  and  Rubesin  (ref.  3).  For 
small  pressure  gradients.  Low  (ref.  4)  has,  by  a  perturbation  analysis, 
treated  the  general  problem  of  the  isothermal  surface.  With  the  intro¬ 
duction  of  pressure  gradients  of  arbitrary  magnitude ,  other  restrictions 
become  necessary.  The  assumption  of  constant  fluid  properties  (density, 
viscosity,  etc.),  for  example,  leads  to  the  greatest  simplification  - 
the  separation  of  the  momentum  and  energy  equations.  With  this  assumption, 
for  a  special  case  of  a  decelerating  stream,  Howarth  (ref.  5)  has  obtained 
a  Beries  solution  to  the  momentum  equation.  The  introduction  of  a  simi¬ 
larity  concept  (that  the  velocity  or  temperature  profiles  may  always  be 
expressed  in  terms  of  a  single  parameter)  leads  to  a  power-law  free-stream 
velocity  distribution.  The  momentum  equation  of  this  problem  was  first 
solved  by  Falkner  and  Skan  (ref.  6),  whose  calculations  were  then  im¬ 
proved  by  Hartree  (ref.  7);  the  energy  equation  was  later  treated  by 
Eckert  (ref.  8)  and  others  (refs.  9  and  10).  For  the  same  problem  the 
restriction  of  constant  fluid  properties  may  be  removed  by  alternatively 
requiring  that  the  Mach  number  be  essentially  zero  (ref.  11)  or  that  the 
Mach  number  and  the  heat  transfer  be  limited  to  small  values  (ref.  12). 

Illingworth  (ref.  13)  and  Stewartson  (ref.  14)  have  demonstrated 
that,  for  an  insulated  surface  in  a  fluid  with  a  Prandtl  number  of  1.0, 
any  compressible  boundary- layer  problem  may  be  transformed  to  a  corre¬ 
sponding  problem  in  an  incompressible  fluid;  the  earlier  solutions  thus 
become  applicable  to  certain  compressible  problems.  For  the  case  of 
heat  flux  across  the  surface,  the  transformation  of  Stewartson  (ref.  14) 
with  the  concept  of  similarity  introduced  leads  to  a  set  of  nonlinear 
ordinary  differential  equations  previously  quoted  (ref.  14),  but  unsolved. 
Solutions  to  this  set  of  equations,  which  are  presented  herein. 


^The  principal  developments  of  this  paper,  which  is  part  of  the  Doc¬ 
toral  Dissertation  of  the  senior  author  (ref.  l),  were  carried  out  under 
the  stimulus  and  guidance  of  Professor  Luigi  Crocco  and  the  sponsorship 
of  the  Daniel  and  Florence  Guggenheim  Foundation.  The  final  analysis 
and  the  computations  were  completed  at  the  NACA  Lewis  laboratory  during' 
the  Spring  of  1954. 
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are  applicable  to  flows  at  arbitrary  Mach  number,  pressure  gradients  of 
arbitrary  magnitude  (but  of  a  form  consistent  with  the  requirements  of 
similarity),  and  arbitrary  but  constant  wall  temperature.^ 

Since  free -stream  velocity  distributions  of  the  form  required  by 
similarity  are  not  generally  encountered  in  practice,  the  utility  of 
these  solutions  is  principally  as  follows:  (l)  the  effects  of  pressure 
gradient,  wall  temperature,  and  Mach  number  may  be  viewed  qualitatively; 
(2)  the  results  may  be  used  as  a  check  on  any  approximate  method  (such 
as  a  Karman-Pohlhausen  method)  for  reliability;  (3)  the  flow  to  be  solved 
may  be  divided  intuitively  into  segments  and  the  solution  for  each 
segment  may  be  matched  by  some  arbitrary  technique;  or  (4)  the  results 
may  be  used  to  construct  a  new  simple  method  (of  the  integral  type)  for 
the  calculation  of  the  laminar  compressible  boundary  layer  with  heat 
transfer.  This  latter  analysis  has  been  carried  out,  utilizing  the 
solutions  herein  given,  and  is  presented  in  reference  1. 


STEWARTSOITS  EQUATIONS 
Boundary -Layer  Equations 

The  equations  of  the  steady  two-dimensional  compressible  laminar 
boundary  layer  for  perfect  fluids  are: 

Continuity: 


Momentum : 


2 

Since  this  writing,  further  calculations,  which  are  closely  re¬ 
lated  to  the  present  investigation,  have  been  published  by  Levy  (ref. 

15).  Solutions  to  the  equations  treated  herein  were  obtained  .in  that 
report.  The' present  investigation  includes  ranges  of  variables  not 
treated  in  ref.  15:  for  example,  favorable  pressure  gradients  applicable 
to  supersonic  nozzles  and  values  of  adverse  pressure  gradients  including 
that  causing  separation.  For  adverse  pressure  gradients,  the  problems  of 
uniqueness  and  multiple  solutions  are  also  considered  in  some  detail. 

The  solutions  of  ref.  15  were  obtained  by  means  of  a  differential  analyzer, 
whereas  the  present  solutions  were  obtained  by  digital  calculation  and 
are  presented  in  tabular  form. 
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Energy: 


dh  dh  Sp  d  At  dh\ 

pu  3x  +  pv  3y  =  c)y(pr  Zy) 

All  symbols  are  defined  In  appendix  A. 

The  viscosity  law  to  be  assumed  is 


/du\2 

+ 


(3) 


iL  =  x  * 
p0  * 


(4) 


0 


Equation  (4)  is  of  the  form  taken  by  Chapman  and  Rubesin  (ref.  3),  ex¬ 
cept  that  the  reference  conditions  (v-q^q)  are  free-stream  stagnation 
values,  since  in  the  presence  of  pressure  gradient  the  local  External" 
values  are  not  constant  along  the  outer  edge  of  the  boundary  layer.  The 
constant  X  is  used  to  match  the  viscosity  with  the  Sutherland  value 
at  a  desired  station.  If  this  station  Is  taken  to  he  the  surface, 
assumed  to  he  at  constant  temperature,  the  result  Is 


+  ksu\ 
+  ksu ) 


(5) 


where  ksu  =  Sutherland's  constant  (for  air,  kBU  =  216°  R) .  The  vis¬ 
cosity  law  of  equations  (4)  and  (5)  was  demonstrated  to  be  adequate  for 
a  flat  plate  (ref.  3)  by  comparison  with  the  more  exact  calculations  of 
reference  2.  In  the  present  case  no  such  comparison  is  available. 


v 


w 


Stewartson's  Transformation 


A  slight  modification  of  Stewartson's  transformation  may  be 
written 

.  ae  Pe  , 
dX  =  X - dx 

a0  p0 


j  -y  P  ae  . 
dY  =  —  r2-  dy 
p0  a0 


> 


(6) 


U  =  \|rY 
Vs  -t. 


J 


c 


t 
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where  the  stream  function  is  defined  by 

7  p0 


The  transformed  quantities  are  now  represented  by  upper-case  letters 
(X,Y,U,V),  and  the  subscript  e  refers  to  local  conditions  at  the  outer 
edge  of  the  boundary  layer  (external).  The  subscript  0  refers  to 
free-stream  stagnation  values.  From  the  preceding  transformation,  a 
useful  relation  between  the  transformed  and  physical  velocities  is 

n  S° 

U  =  —  u. 

ae 

If  equations  (4)  and  (6)  are  applied  to  the  boundary- layer  equations 
(l).»  (2),  and  (3),  and  if  Pr  and  Cp  are  taken  to  be  constant  (but 
it  is  not  yet  required  that  Pr  =  l),  there  result 

Ux  +  VT  =  0  (7) 


UUX  +  WY  —  UqUq  (l  +  S)  +  vqUyy 


tto  .  ™  ..  )  YY  1  -  Pr 

USX  +  VSY  =  vQ(  - 


-  im2  \r  " 
2  Me  \  /V\2 

7-1  m2  J  l  U  J 

2  Me/  '  ✓  YY 


where  the  enthalpy  function  S  is  defined  for  convenience  as 


S  =  >T  "  1 
h0 


and  h8  is  the  local  stagnation  enthalpy. 
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The  boundary  conditions  applicable  to  the  system  (7)  to  (9) 

are : 


•N 

u(x,o)  =  o 

V(X,0)  =  0 

S(x,0)  =  Sw  or  jj|  (X,0)  =  y 

lim  S  =  0 
Y-+00 

lim  U  =  Ue(X) 

Y-+“> 


(ID 


The  solution  S  =  0  and  the  resultant  continuity  and  momentum  equations 
(7)  and  (8)  make  up  the  extremely  useful  correlation  developed  by 
Stewartson  between  compressible  and  incompressible  boundary  layers  on 
insulated  surfaces  with  Pr  =  1.  Another  special  case  is  that  of 

XJ0^.  =  0.  Then,  if  Pr  =  1,  the  relation  S  =  Sw^l  -  satisfies 

equation  (9);  this  is  Crocco's  integral  of  the  energy  equation  for  the 
flat  plate  (ref.  2). 


v 


V 


Similarity  Requirements 

When  a  pressure  gradient  exists  and  the  surface  is  not  insulated, 
it  is  necessary  to  find  a  means  of  solving  the  system  (7)  to  (9)  subject 
to  the  boundary  conditions  (ll).  To  this  end,  the  question  will  be 
asked:  Under  what  conditions  can  this  system  be  reduced  to  a  system  of 
ordinary  differential  equations  by  the  assumption  that  the  boundary- 
layer  profiles  are  functions  of  a  similarity  variable  t)  and  that  the 
wall  temperature  is  constant?  This  question  may  be  resolved  by  inserting 
the  following  assumed  relations  into  the  system  (7)  to  (9)  and  observing 
the  conditions  required  for  obtaining  ordinary  differential  equations: 

*  =  AxV  f^ 

Y  =  BX^It) 

S  =  S(tj) 

where  A,B,a,b,p,  and  q  are  undetermined  constants.  This  procedure 
has  been  carried  out  by  Li  and  Nagamatsu  (ref.  16)  for  Pr  =  1.  In 
that  analysis  it  was  concluded  that  four  classes  of  similar  solutions 


(12) 
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are  possible.  It  has  been,  pointed  out  (ref.  17)  that  three  of  these 
four  classes  can  be  reduced  identically  to  the  case  requiring  that 


Ue  =  CX® 

while  the  remaining  case  requires  that 


UQ  =  G]_  exp  C2X 


(13) 


(14) 


When  equations  (12)  are  used  in  the  form 

* 

t-  f(n) 


(15) 


IP  + 

1  ue 

II  2 

V 

the  system  of  ordinary  differential  equations  corresponding  to  the 
power-law  velocity  distribution  of  equation  (13)  may  be  written 

f +  ff"  =  p(f ,2  -  1  -  S) 

(7  "  1)M| 


S”  +  PrfS '  =  (1  -  Pr) 


1  +  2-f1^ 


(ff'"  +  f"2) 


2m 


(16) 


and  the 


The  pres sure -grad lent  parameter  P  is  defined  as  p  = 

m  +  1' 

velocity  ratio  Is  U/UQ  =  u/u0  =  f ' ,  where  primes  denote  differentiation 
with  respect  to  t). 


The  houndary  conditions  are: 


f  (0)  =  f  '(0)  =  <0 
s(0)  =  S* 

lim  f '  =  1  > 

*1-** 


lim  S  =  0 

J 


(17) 


Since  Mq  may,  in  general,  he  a  function  of  x;  the  right  member  of 
the  energy  equation  is  not  yet  dimensionally  consistent  with  the  left 
member  for  arbitrary  Me  and  Pr. 
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It  was  showfr  in  reference  17  that,  for  the  exponential  case 
(eq.  (14))  with  Cg  >  0,  the  system  (7)  to  (9)  can  he  reduced  to  the 
ordinary  differential  equations  (16),  hut  with  0  =  2.  For  C£<0,  the 
f™  term  in  equations  (16)  is  replaced  hy  -f ,H .  In  this  case,  with 
S  =  0,  it  can  he  shown  that,  because  of  the  sign  of  the  f m  term,  no 
solution  is  possible  in  which  the  velocity  ratio  approaches  its  boundary 
condition  smoothly.  A  question  is  thus  raised  as  to  the  validity  of 
any  possible  solution  for  Cg  <  0  regardless  of  the  value  of  S.  For 
the  remainder  of  this  paper  this  class  will  be  omitted  from 
consideration. 


Corresponding  analyses  for  incompressible  flow,  including  condi¬ 
tions  for  similarity  and  the  case  of  the  exponential  free-stream  veloc¬ 
ity,  have  been  made  by  Mangier  (ref.  18)  and  Goldstein  (ref.  19), 
respectively.  As  previously  mentioned,  the  right  member  of  the  energy 
equation  (16)  must  be  zero  or  a  function  of  tj  to  be  consistent  with 
the  left  member.  This  may  be  achieved  in  the  following  ways:  (l)  the 
external  Mach  number  may  be  a  constant  other  than  zero,  (2)  the  external 
Mach  number  may  be  zero,  (3)  the  Prandtl  number  may  equal  1,  (4)  the 
(7  -  i)m§  I 

2  corresponding  to  hypersonic  flow,  or  (5) 


factor 


1  + 


2  e 

the  ratio  of  specific  heats 


7  may  equal  1. 


The  case  of  constant  external  Mach  number  is  the  flat-plate  problem 
(0  =  0)  and,  the  solution  to  the  momentum  equation  being  known,  the 
energy  could  be  integrated  directly.  The  flat-plate  problem  has  already 
been  solved  with  great  accuracy  and  completeness  by  Crocco  (ref.  2). 

If  the  pressure  gradient  is  small  enough,  it  may  be  reasonable  to  con¬ 
sider  Mq  constant  in  the  energy  equation  in  spite  of  the  gradient, 
but  to  retain  the  pressure-gradient  parameter  in  the  momentum  equation. 
However,  this  problem  is  treated  more  completely  by  the  analysis  of 
reference  4. 


The  case  Mq  =*  0  (with  arbitrary  0)  produces  the  equations  of 
Levy  and  Seban  (ref.  20).  In  that  analysis  approximate  solutions  were 
obtained  by  the  assumption  of  simple  forms  for  the  velocity  and  tempera¬ 
ture  profiles  which  contained  undetermined  coefficients.  These  coeffi¬ 
cients  were  then  evaluated  by  use  of  the  boundary  conditions.  Because 
the  actual  profiles  cannot  be  simply  represented,  this  method  is  not 
reliable  in  some  ranges  even  if  the  Mach  number  is  nearly  zero.  Brown 
and  Donoughe  (ref.  ll)  also  considered  the  low  Mach  number  problem  with 
variable  fluid  properties  and  Prw  =  0.7.  The  system  of  equations 
encountered  in  that  analysis  is  much  more  complicated  than  the  present 
system  because  of  the  power-law  viscosity,  conductivity,  and  specific- 
heat  relations  used.  These  refinements  do  not  alter  the  effects  of 
omitting  the  viscous-dissipation  and  compressive -work  terms,  which  may 
be  significant  at  higher  Mach  numbers. 


v 
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The  case  of  hypersonic  flow  requires  the  introduction  of  the  effects 
of  displacement  thickness  upon  pressure  gradient,  such  as  have  been 
evaluated  by  Lees  and  Probstein  (ref.  21),  for  example.  This  case  will 
not  be  treated  herein. 

The  possibility  of  assuming  7  =  1  does  not  simplify  the  equations 
beyond  the  assumption  of  Mach  number  zero.  For  most  gases,  the  assump¬ 
tion  of  7=1  is  physically  unreasonable.  Therefore,  this  case  does 
not  appear  to  warrant  further  consideration. 

If  strong  pressure  gradients  and  reasonably  high  Mach  numbers  are 
to  be  considered,  it  thus  appears  desirable  to  restrict  the  similarity 
system  to  Pr  =  1,  with  the  result  that 

f +  ff"  =  p(f ,2  -  1  -  S)  (18a) 

S"  +  fS'  =0  (18b) 

with  the  boundary  conditions  (17).  Equations  (18)  were  derived  by 
Stewartson  by  assigning  similarity  relations  corresponding  to  (15)  to 
the  system  (7)  to  (9)  with  Pr  =  1;  however,  no  solution  was  indicated. 

The  comparison  between  assuming  that  Mq  =  0  (case  (2))  or  that 
Pr  =  1  (case  (3))  may  perhaps  be  indicated  by  examination  of  the 
solutions  to  the  insulated  flat-plate  problem,  which  include  effects 
of  both  Prandtl  number  and  Mach  number  (ref.  2).  If  Mq  =  0,  the 
viscous-dissipation  and  compressive-work  terms  are  omitted  in  equation 
(3).  Then  the  predicted  temperature  profile  is  a  constant,  rather  than 
the  correct  variation  from  free-stream  static  to  recovery  temperature 
at  the  wall.  However,  if  Pr  »  1  is  assumed,  a  constant  stagnation 
temperature  is  predicted,  rather  than  the  actual  slight  variation  in  this 
quantity.  The  latter  discrepancy  is  small  compared  with  the  former. 


METHOD  OF  SOLUTION 

Equations  (18)  with  boundary  conditions  (17)  comprise  the  system 
to  be  solved  for  the  dependent  variables  f(q)  and  S(tj).  Because  of 
the  nonlinearity  of  the  system,  its  high  order  (fifth),  and  its  classi¬ 
fication  as  a  "two-point  boundary- value  problem,"  no  standard  integration 
methods  will  yield  results  expressible  in  closed  form.  Methods  appli¬ 
cable  to  equations  of  this  type  may  be  classified  as  either  (l)  forward 
integrations  or  (2)  integrations  by  methods  of  successive  approximations. 

By  "forward  integration"  is  meant  the  progressive  integration  of 
the  equations  from  one  (initial)  boundary  to  the  other.  For  this 
purpose  several  sets  of  initial  values  of  the  derivatives  are  assumed. 
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Then  the  final  boundary  values  obtained  are  compared  with  those  speci¬ 
fied  and,  after  interpolation  of  the  initial  values,  this  trial-and- 
error  process  is  repeated  until  the  final  boundary  conditions  are  satis¬ 
fied.  The  integrations  may  be  carried  out  by  the  use  of  either  an 
analog  computer  (mechanical  or  electrical)  giving  continuous  integrals 
or  by  digital  computations  involving  finite-difference  integration. 
Although  generally  applicable,  a  disadvantage  associated  with  forward 
integration  of  nonlinear  equations  is  the  lack  of  any  inherent  conver¬ 
gence  mechanism.  Thus,  the  approach  to  the  correct  initial  values 
depends  almost  entirely  on  the  intuition  and  experience  of  the  one 
performing  the  calculations.  This  method  is  particularly  troublesome 
for  a  problem  with  more  than  one  dependent  variable  since  evidence  for 
the  fitness  of  a  given  initial  value  may  be  obscured  by  a  poor  selection 
of  the  corresponding  initial  value  of  another  dependent  variable. 
Furthermore,  when  an  analog  computer  is  employed  the  accuracy  is  general¬ 
ly  limited,  particularly  for  nonlinear  equations  where  in  certain  regions 
the  results  tend  to  be  highly  sensitive  to  the  chosen  initial  values. 

If  digital  computation  is  utilized  to  obtain  a  desired  degree  of  accuracy, 
the  procedure  may  become  excessively  tedious. 

Successive  approximation  methods  generally  assume  an  entire 
function  for  the  dependent  variables  (satisfying  as  many  of  the  boundary 
conditions  as  possible)  rather  than  only  the  initial  derivatives.  Then, 
by  use  of  the  differential  equations,  a  procedure  is  developed  for 
estimating  the  error  as  a  function  of  the  independent  variable(s).  This 
error  is  applied  to  the  original  choice  and  the  process  is  repeated 
until  satisfactory  convergence  occurs.  An  example  of  a  method  of 
successive  approximation  is  Picard's  method. 

A  difficulty  shared  by  both  these  methods  arises  when  the  range  of 
integration  is  infinite.  Then  it  is  necessary  to  decide  upon  a  finite 
value  of  the  independent  variable  at  which  the  boundary  conditions  may 
be  approximately  satisfied  and  the  degree  to  which  they  may  be  satisfied. 
This  suggests  the  desirability  of  changing  to  an  independent  variable 
so  that  only  a  finite  range  of  integration  is  required.  In  the  present 
problem  this  change  of  variables  can  be  achieved  by  following  a  method 
used  by  Crocco  for  the  solution  of  the  compressible  flat-plate  boundary 
layer  (ref.  2).  The  concept  is  advanced  that  the  velocity  is  a  more 
suitable  independent  variable  Bince  it  is  bounded.  This  concept  leads 
to  a  set  of  equations  conveniently  handled  by  a  method  of  successive 
approximations . 


Transformation  to  Velocity  Plane 


To  accomplish  the  transformation  to  the  velocity  ratio  f'  as  the 
independent  variable,  the  following  identity  may  be  used: 


d  _  d 
drj  —  1  df ' 


(19) 


» 
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This  Identity  may  be  applied  to  f "  and  f 

f =  f”  ” 
df ' 

°0  v o  f 


as  follows: 

■N 


(20) 


where  the  dummy  variable  of  integration  is  and  f"(£)  represents 
the  functional  relationship  between  f"  and  f ' ,  that  is,  f "(f). 
The  primes  continue  to  denote  differentiation  with  respect  to  tj. 


in 


Inserting  equations  (20)  into  the  momentum  equation  (18a)  results 


df” 
df ' 


+  P 


(21) 


Hi 


which  satisfies  the  following  condition  at  f f  =  0  required  b y  the 
momentum  equation: 


f^  =  -  p(l  +  Sy)  (22) 

Now,  if  equation  (21)  is  integrated  once  with  respect  to  f f  and  if 
the  limits  of  integration  are  chosen  so  that  (f")^^  =  0,  the  result 
is 


By  inverting  the  order  of  integration  (or  "by  integrating  by  parts) 
the  double  integral  may  be  reduced  to  two  single  integrals ,  resulting 
in: 


fj+i(f')  = 


r 1  iizzim  -  ft  f 1  s2- i-s^c) 
Jf  fjttr  p  1,  -i  id) 


d? 


(24) 


Equation  (24)  is  the  form  of  the  momentum  equation  as  it  will  be  used 
in  this  report.  The  subscript  j  is  the  iteration  number  in  the 
method  of  successive  approximations. 


A  corresponding  form  of  the  energy  equation  is  obtained  by  writing 
equation  (18b)  as 
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and  integrating  with  respect  to  tj,  to  get 

In  S'  =  -  ^  f  drj  +  constant 
Equation  (18a)  may  "be  written 

(f2  -  1  -~s) 


(25) 


-P  I  Tt 

f  dTl  =  -  p-  dTl  +  0 


f" 


dr) 


_  _  df"  p  (f’2  -  1  -  S) 


f" 


(f'T 


df 


Substitution  of  this  expression  into  equation  (25)  results  in 


In  S’ 


g2  -  1  -  S(g) 
[f"(g)  ]2 


dg  +  constant 


or  the  equivalent  expression 


where 


S'  =  -  C3f"j(f) 


J(g)  =  exp 


-0 


g*  -  i  -  s(eL) 

[f-Ui)]2 


(26) 


If  this  expression  is  integrated  once  again  and  the  boundary  con¬ 
ditions  S(0)  -  (S)^  t_i  =  0  are  required,  the  result  is 


(27) 


Inspection  of  equations  (24)  and  (27)  indicates  that  the  integrals 
to  he  evaluated  are  singular,  or  indeterminate,  at  the  upper  limit.  To 
evaluate  these  integrals,  closed-form  expressions  must  he  obtained  for 
the  integrands  in  this  range.  This  requires  knowledge  of  the  solution 
of  the  system  (18)  for  large  r\  (near  f’  =  l) .  This  "asymptotic  solu¬ 
tion"  and  its  development  are  given  in  appendix  B.  The  results  show 
that  equation  (24)  can  he  used  in  its  present  form,  hut  that  equation 
(27)  must  he  modified  to 


eJj(l-e)  + 


ejj(l-e)  + 


1-6 


! 


Jj(g)dg 


Jj(g)ag 


v 


(28) 


v 
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where  e  is  an  arbitrary  small  quantity  (e  «  1).  In  this  form  the 
singularity  has  been  removed.  Equations  (24)  and  (28)  comprise  the 
system  used  in  the  present  investigation.  The  convergence  of  this  sys¬ 
tem  is  discussed  in  appendix  C,  and  the  method  of  calculation  in  appen¬ 
dix  D. 


PROPERTIES  OF  SOLUTIONS 

In  the  following  sections  the  Solutions  obtained  in  this  study  are 
presented  and  their  properties  are  discussed.  The  two  parameters  de¬ 
fining  a  case  are  S-y  and  0.  The  enthalpy  function  evaluated  at  the 
wall  Sy-  determines  the  wall  temperature  through  the  relation 

tw  =  t0(l  +  S*)  (29) 

Thus ,  Sw  =  -1  corresponds  to  a  wall  temperature  of  absolute  zero, 
and  3^=1  corresponds  to  a  wall  at  twice  the  free-etream  stagnation 
temperature.  The  case  Sy-  =  0  corresponds  to  a  wall  at  the  free-stream 
stagnation  temperature,  which  for  Pr  =  1  is  the  case  of  an  insulated 
surface . 

The  pressure-gradient  parameter  p  is  related  to  the  exponent  m 
of  the  velocity  distribution  in  the  transformed  plane  Ue  =  CX™  through 
the  relation 


For  a  velocity  distribution  of  this  form,  m  can  be  represented  as 

°  ■  GrM  S0  31  (30) 

It  is  apparent  that  p  <  0  (m  <  0)  corresponds  to  an  unfavorable 
gradient;  p  =  0  (m  =  0)  corresponds  to  flat-plate  flow;  and  p  =  2 
(m  =  »)  corresponds  to  an  infinitely  favorable  pressure  gradient. 
Stewartson  (ref.  14)  has  shown  that  p  =  1  (m  =  l)  corresponds  to  flow 
in  the  immediate  vicinity  of  a  stagnation  point  for  two-dimensional  flow, 
as  in  the  incompressible  case.  It  can  be  shown  that  the  case  of  a  stag¬ 
nation  point  in  axisymmetric  flow  can  be  transformed  to  the  solution 
for  p  =  l/2  (ref.  22).  An  approximate  method  for  relating  p  to  more 
general  physical  flows  is  given  in  reference  1.  Values  of  p  of  the 
order  of  magnitude  ±0.3  correspond  to  flows  over  supersonic  wings,  and 
a  typical  nozzle  with  an  exit  Mach  number  of  about  2.5  might  produce  a 
value  of  P  of  about  1.5.  In  the  present  investigation,  solutions  are 
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found  for  pressure  gradients  ranging  from  that  causing  separation  to  the 
infinitely  favorable  gradient  and  for  wall  temperatures  from  absolute 
zero  to  twice  free-stream  stagnation  temperature.3 

All  solutions  are  presented  in  tabular  and  graphic  form.  Table  I 
shows  the  values  of  f,  f',  f ",  S,  and  S'  tabulated  against  r\.  From 
these  values  and  equations  (18)  the  quantities  f"'  and  S"  can  be  easily 
calculated.  Table  II  presentb  a  summary  of  the  values  of  fy  (related 
to  wall  shear)  and  'S^  (related  to  heat  transfer)  from  table  I,  as  well 
as  the  Reynolds  analogy  parameter  CfRew/Ru,  which  represents  the  ratio 
of  shin-friction  to  heat-transfer  effects.  Certain  other  quantities  of 
interest  cannot  be  tabulated  in  general,  but  can  be  easily  calculated 
from  the  following  formulas : 


Static -temperature  ratio: 


X^-iM^(l  +  S)  -  2L=_iM§  f 


2 


(31) 


or,  with  the  static  temperature  t  referred  to  the  free-stream  stagna¬ 
tion  temperature  tQ, 


Flux  density: 


(32) 


Pu 

peue 


f ' 

2-j-i  Mijd  +  S) 


7  -  1  M2f ,2 
2  6 


(33) 


Uniqueness 


For  p  <  0,  Sy  =  0,  Hartree  (ref.  7)  first  observed  that  the  bound¬ 
ary  conditions  (17)  are  not  sufficient  to  determine  a  unique  solution. 


3lt  should  be  noted  that  all  but  one  of  the  presented  solutions  for 
Sy  =  0  are  thuse  first  obtained  by  Hartree  (ref.  7)  for  the  problem  of 
Falkner  and  Shan  (ref.  6).  As  a  further  chech  on  the  present  method, 
the  solutions  for  P  =  1.6  and  2.0  with  Sy  =  0  were  obtained  independ¬ 
ently  in  the  present  investigation;  these  values  agree  very  well  with 
those  of  Hartree. 


t 
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Thus,  there  is  not  a  unique  value  of  for  a  given  p.  In  studying 

the  uniqueness,  it  is  useful  to  consider  the  following  expression  for 
velocity  ratio  (for  any  S^)  valid  for  large  tj: 

f  ’  =  1  +  ^(ti  -  x)-(20+l)  +  —■  (t]  -  x)~ijexp£  1^  ■"  J  +  02  ( T)  -  x)2(3 

(34) 

where  Oq_,  a g,  a^,  and  x  are  integration  -constants  (see  appendix  B). 

Hi  case  of  S^-  =  0,  is  also  equal  to  zero;  however,  this  does  not 

change  the  uniqueness  problem,  which  is  independent  of  wall  temperature. 
For  p  >  0,  ag  is  necessarily  zero  in  order  to  satisfy  the  boundary 
condition  lim  f f  =  1,  For  continuity  in  .p,  Hartree  then  selected  the 

T)-*» 

asymptotic  solution  with  cxg  =  0  for  p  <  0. 


Another  important  result  of  the  asymptotic  solution  is  ‘that  the 
integral 

X”  (x ' pu )  dr),  related  to  the  displacement  thickness, 

can  he  shown  to  become  infinite  for  a^j  ^  0.  This  result  is  contrary  to 
the  concept  of  a  thin  layer  outside  of  which  the  viscous  effects  may  he 
neglected.  A  further  effect  of  the  otp  term  on  the  solution  can  he 
observed  by  examination  of  the  dimensionless  quantity  f '"/ff"  (sug¬ 
gested  by  Professors  L.  Crocco  and  L.  Lees),  in  which  f"  represents 
the  net  viscous  forces  acting  on  the  fluid  element  and  f"  is  propor¬ 
tional  to  the  velocity  gradient  (shearing  flow).  It  can  be  shown  that 
for  ag  =  0 


while  for  cig  ^  0 


Solutions  with  ctg  =  0  retain  the  numerator  and  the  denominator  of  the 
ratio  -f /ff "  to  the  same  order  of  magnitude,  while  if  Og  is  dif¬ 
ferent  from  zero  a  solution  results  wherein  the  magnitude  of  the  net 
viscous  forces  in  the  asymptotic  region  is  small  compared  with  the  mag¬ 
nitude  of  the  shearing  flow  set  up  by  their  action.  Thus,  in  order  to 
retain  both  effects  of  viscosity  to  the  same  order  of  magnitude,  ctg 
must  be  taken  equal  to  zero,  as  was  done  by  Hartree.  The  solution  thus 
obtained  will  be  termed  the  "viscid"  solution. 
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Another  feature  of  solutions  with  different  from  zero  is  the 

analytical  result  that  the  velocity  ratio  in  the  outer  portion  of  the 
"boundary  layer  may  exceed  unity.  For  example,  if  is  not  zero,  equa¬ 
tion  (34)  shows  that  for  large  t)  the  term  of  the  velocity  ratio 

expression  is  dominant,  and  thus  (f f  -  l)  is  necessarily  of  the  same 
sign  as  a^.  That  is,  for  positive  a 2,  the  velocity  ratio  approaches 
unity  from  above;  this  phenomenon  will  be  termed  "velocity  overshoot." 
Since,  for  a  given  p  and  Sw  in  this  range,  each  of  these  various 
solutions  has  associated  with  it  a  different  set  of  values  of  f "  and 

Vf 

S^,  one  of  these  parameters,  say  f^,  can  be  conveniently  used  in  place 
of  02  to  identify  the  various  solutions.  This  infinite  set  of  solutions 
can  be  represented  as  in  sketch  (a)  for  a  typical  (cold  wall)  case. 


It  is  seen  that  there  are  a  maximum  and  a  minimum  shear  (represented  by 
f^)  and  heat  transfer  (represented  by  S^J.)  that  can  satisfy  the  equations 
without  incurring  velocity  overshoot.  These  distinct  solutions  (circled 
points  in  sketch  (a))  correspond  to  og  =  o,4  the  viscid  solutions;  that 


4In  the  evaluation  of  the  singularities  of  the  integrals  required 
for  the  method  of  successive  approximations,  Og  was  taken  to  be  zero. 
Hence,  solutions  for  ocg  ^  0  were  obtained  by  forward  integration 
(appendix  D),  although  the  numerical  values  of  ctg  were  not  determined. 
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with  the  lower  shear  Is  designated  the  "lower-hranch "  solution.  The 
behavior  of  the  calculated  family  of  solutions  is  presented  in  figure  1 
for  Sy.  =  -0.8  and  3  =  -0.325,  -0.3285,  and  -0.336.  For  a  given  value 
of  Sy-,  as  3  is  decreased  the  two  viscid  solutions  approach  each  other. 
At  a  value  of  3  to  he  designated  3^^,  these  two  solutions  become 
identical  and  for  3  <  no  viscid  solution  exists.  For  negative  3, 

only  the  viscid  solutions  will  be  considered  in  the  remainder  of  this 
report . 

With  regard  to  the  physical  significance  of  the  double  solution,  it 
may  be  noted  that  for  adverse  pressure  gradients  (3  <  0)  a  real  flow 
cannot  completely  reproduce  the  similar  solution  because  Ue(0)  =  °°  would 
be  involved.  However,  a  pressure  field  can,  in  principle,  be  applied  to 
a  developing  boundary  layer  so  that,  after  a  phase  of  adjustment,  the 
boundary  layer  would  approach  one  of  the  similar  solutions  with  3  <  0 
and  stay  quite  close  to  it  thereafter.  It  seems  reasonable  to  believe 
that,  depending  on  the  way  the  pressure  field  is  applied,  one  solution 
or  the  other  corresponding  to  the  same  3  could  be  approached  after 
different  adjustment  phases.  This  result  is  exactly  what  Clauser  (ref. 

23)  has  found  in  his  experimental  work  on  Bimilar  turbulent  boundary- 
layer  flows. 


Telocity  and  Temperature  Profiles 

The  velocity  and  enthalpy- function  profiles  obtained  from  the  tabu¬ 
lated  solutions  are  presented  as  functions  of  tj  in  figures  2  and  3, 
respectively.  The  distance  y  normal  to  the  surface  in  the  physical 
plane  is  related  to  the  similarity  variable  tj  through  equations  (6) 
and  (15),  and  may  be  expressed  as 


where  t/tQ  is  given  by  equation  (32). 

Telocity  overshoot.  -  The  velocity  profiles  shown  in  figure  2  indi¬ 
cate  that  for  a  given  wall  temperature  the  initial  slope  decreases  as 
the  pressure  gradient  becomes  less  favorable.  For  adverse  pressure 
gradients  an  inflection  point  occurs  within  the  boundary  layer -and  moves 
outward  as  the  gradient  becomes  more  adverse.  The  velocity  ratio  varies 
monotonically  from  zero  to  the  final  value  of  1.0  except  for  the  cases  of 
favorable  pressure  gradients  with  heated  walls .  Then  the  velocity 
ratio  in  the  outer  portion  of  the  boundary  layer  reaches  a  maximum  value 
greater  than  1.0  before  returning  to  its  final  value  of  1.0  This  type 
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of  velocity  overshoot  was  also  obtained  in  the  investigation  of  reference 
11  for  favorable  pressure  gradients  with  heated  walls  and  is  to  be  dis¬ 
tinguished  from  that  associated  with  the  nonunique  inviscid  solutions 
which  occur  only  for  adverse  pressure  gradients.  When  the  wall  is 
heated  in  a  favorable  pres sure -grad lent  flow,  the  density  within  certain 
layers  of  the  boundary  layer  is  lowered  so  that,  in  spite  of  the  viscous 
retardation,  the  flow  is  accelerated  more  than  the  external  flow  by  the 
external  pressure  forces.  Thus,  a  velocity  greater  than  the  external 
velocity  may  be  obtained. 

This  phenomenon  can  be  established  by  examination  of  equation  (34) 
and  the  corresponding  asymptotic  expression  for  the  enthalpy  function 
(appendix  B) : 


S  =  ajOl 


2 


(36) 


For  favorable  pressure  gradients,  =  0  as  previously  mentioned.  Then, 
the  cxg  term  in  equation  (34)  is  dominant  for  large  t| .  Thus,  (f '  -  1) 
and  oc^  are  of  the  same  sign.  Hence,  for  a  heated  wall  (a^  positive, 
eq.  (36)),  the  velocity  ratio  must  approach  1.0  from  above. 


Stagnation-temperature  profiles.  -  Figure  3  shows  that  for  Pr  =  1, 
the  stagnation  temperature  varies  monotonically  across  the  boundary 
layer  from  the  wall  value  to  the  free-stream  value.  For  favorable  pres¬ 
sure  gradients  with  a  cold  wall,  there  is  small  variation  with  P  of 
this  distribution.  The  variation  becomes  more  pronounced  with  an  in¬ 
crease  in  wall  temperature. 

Boundary- layer  thickness.  -  The  velocity  profiles  (fig.  2)  indicate 
that  the  boundary  layer  thickens  as  the  wall  shear  stress  diminishes. 
Also,  for  a  given  value  of  the  pres sure -grad lent  parameter  p,  the  boun¬ 
dary  layer,  when  considered  in  terms  of  tj,  thickens  as  the  wall  tempera¬ 
ture  is  lowered.  However,  in  the  physical  plane  (in  terms  of  y)  be¬ 
cause  of  the  relation  between  y  and  t)  (eq.  (35))  the  trend  is  Just 
the  opposite.  This  emphasizes  the  necessity  for  careful  consideration 
of  the  relation  between  the  transformed  quantities  and  their  physical 
counterparts. 


The  thermal  boundary  layer  also  thickens  as  separation  is  approached. 
The  relative  thicknesses  of  the  dynamic  and  thermal  boundary  layers  may 
be  conveniently  observed  from  a  plot  of  S  against  f'  (fig.-  4).  Then 
if  a  fixed  fraction  of  S^,  say  0.99,  is  chosen  to  define  the  thermal- 
layer  thickness  and  if  the  same  value  of  velocity  ratio  is  taken  to  de¬ 
fine  the  dynamic  layer,  it  can  be  seen  that,  regardless  of  wall  tempera¬ 
ture,  the  thermal  layer  is  thicker  than  the  dynamic  layer  for  favorable 
gradients  and  thinner  for  adverse  gradients. 


w 
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For  Pr  <  1  the  relative  magnitude  of  the  dynamic  thickness  to 
the  thermal  thickness  will  he  decreased,  since  the  Prandtl  number  re¬ 
presents  the  ratio  of  viscous  to  thermal  effects  in  the  fluid. 


Shear  and  Skin  Friction 

The  shear  distribution  in  the  boundary  layer  is  presented  in 
figure  5,  where  f"  is  plotted  as  a  function  of  r\,  The  shear  function 
f"  is  related  to  the  shear  stress  t  through  the  expression 


27-1 

- 
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\Jm  +  1 

V  2  vox 

(37) 


For  (3  >  0  the  maximum  shear  is  at  the  wall,  whereas  for  (3  <  0  the 
point  of  maximum  shear  moves  increasingly  outward  as  the  pressure 
gradient  becomes  more  adverse. 


The  quantity  that  is  of  primary  interest  in  boundary- layer  calcula¬ 
tions  is  the  shear  stress  at  the  wall  xv,  which  can  be  made  dimension¬ 
less  through  the  definition  of  a  local  skin-friction  coefficient, 
producing  the  relation 


(38) 


The  factor  (l  +  S^)  appears  in  equation  (38)  because  of  the  use  of  p 
in  the  definition  of  Cf.  Although  this  factor  can  be  easily  avoided, 

...  uax 

it  is  used  later  in  evaluating  a  Reynolds  number  Rew  =  — —  suitable 

w 

for  use  in  determining  the  heat  transfer.  An  alternate  form  for 
equation  (38)  is 


=  r a  /m T  1  d  liTx 

2  w  y  2  d  In  x 


(38a) 


It  should  be  noted  that,  in  equation  (38a),  fluid  properties  are  evalu¬ 
ated  at  the  wall  temperature.  If  the  skin-friction  coefficient  and 
the  Reynolds  number  were  to  be  based  on  local  free -stream  fluid 
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properties,  rather  than  on  wall  values,  a  factor  of 


V  ^e  tw 


—  would 


appear  in  the  right  member  of  equation  (38a)  .  When  this  factor  is 

,1/4 

r  vtr  * 

evaluated  using  Sutherland's  viscosity  law,  it  varies  from 


V 


-1/4 


Aw 

\  ® 


to 


— J  ,  depending  on  the  temperature  involved 
te  1 


The  quantity  f^  is  presented  as  a  function  of  (3  and  Sw  in 
figure  6.  It  can  he  seen  that  heating  the  surface  increases  the  sensi¬ 
tivity  of  the  wall  shear  to  pressure  gradient,  while  cooling  the  wall 
has  the  opposite  effect.  A  suggested  physical  interpretation  for  this 
trend  is  related  to  the  effect  of  wall  temperature  on  the  mean  density 
of  the  fluid  within  the  boundary  layer.  For  the  heated  wall,  the 
houndary-layer  density  is  less  than  the  free -stream  density,  rendering 
the  boundary- layer  fluid  more  susceptible  to  free-stream  acceleration 
forces  than  for  the  cold  wall.  Figure  6  shows  further  that  a  linear 
extension  of  the  slope  of  the  curve,  f^  against  P  from  p  =  0  to 
large  positive  p  would  grossly  overemphasize  the  effects  of  favorable 
pressure  gradient;  while  the  same  linear  extension  toward  negative  p 
would  underemphasize  the  effects  of  adverse  gradient. 


In  figure  6(b),  the  two  viscid  solutions,  which  occur  for  adverse 
pressure  gradients  for  a  given  P  and  Sy,  are  plotted.  It  is  seen 
that  two  solutions  are  given  for  even  the  insulated  surface  (Sv  ±  0), 
although  Hartree  reported  only  one.  In  this  case  the  lower-branch 
solution  corresponds  to  negative  wall  shear  stress  (separated  flow), 
which  was  not  considered  in  reference  7.  For  heated  walls  >  0) 
both  solutions  may  be  separated  near  Pmin>  while  for  cooled  walls 
both  solutions  may  be  unseparated  in  this  region.  The  physical  inter¬ 
pretation  of  these  double  solutions  has  been  discussed  in  the  section 
UNIQUENESS . 


* 


Heat  Transfer 


The  variation  of  heat  transfer  across  the  boundary  layer  is  plot¬ 
ted  in  figure  7  in  terms  of  the  derivative  of  the  enthalpy  function 

S’  =  — .  This  quantity  is  related  to  the  stagnation  enthalpy  deriva- 

dr) 

tive  in  the  physical  plane  by  the  expression 
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These  carves  again  indicate  the  thickening  of  the  thermal  layer  as 
separation  is  approached.  Furthermore,  as  separation  is  neared,  the 
zone  adjacent  to  the  surface  where  S’  is  essentially  constant  spreads 
rapidly.  This  is  a  zone  where  the  heat  transfer  is  primarily  by  con¬ 
duction  because  of  the  near  zero  velocities  in  the  neighborhood  of 
the  surface. 

The  values  of  Sf  at  the  surface  (S^)  are  shown  plotted  as  a 
function  of  pressure-gradient  parameter  (3  in  figure  8  for  constant 
wall  temperatures.  Two  facts  are  noteworthy:  (l)  In  the  region  of 
favorable  pressure  gradient,  is  nearly  constant;  (2)  the  heat 
transfer  varies  sharply  near  separation.  From  these  facts  the  addi¬ 
tional  conclusion  may  be  drawn  that,  if  a  linear  extension  of  these 
curves  is  made  with  the  slope  at  (3  =  0,  the  result  will  seriously 
overemphasize  the  effects  of  a  favorable  pressure  gradient  or  heat 
transfer  and  underestimate  the  effects  for  adverse  pressure  gradients. 

A  similar  influence  of  pressure  gradient  on  skin  friction  has  already 
been  noted.  A  comparison  of  figures  6  and  8  indicates  that  the  effect 
of  pressure  gradient  on  heat  transfer  is  smaller  than  the  corresponding 
effect  upon  wall  shear. 

As  with  the  skin  friction,  it  is  convenient  to  define  a  dimension¬ 
less  number  from  which  the  heat  transfer  may  be  determined.  The  Nusselt 
number  is 


m  +  1  d  In  X 
2  d  In  x 


(40) 


The  quantity  (-S^/S  )  is  plotted  in  figure  9  for  constant  wall  tempera¬ 
tures  as  a  function  of  the  pres  sure -grad  lent  parameter  (3.  The  Reynolds 
number  Rew  is  again  defined  in  terms  of  wall  properties. 

Reynolds  analogy.  -  From  expressions  (38a)  and  (40),  a  simple 
modifed  Reynolds  analogy  parameter  is  evaluated  by 


Nu 


(41) 


This  quantity  is  the  reciprocal  of  the  usual  Reynolds  analogy  quantity 
in  order  to  avoid  infinite  values  as  separation  is  approached.  It 
is  plotted  in  figure  10  as  a  function  of  the  pres sure -gradient  parameter 
(3 .  These  curves  resemble  the  curves  (fig.  6)  because  of  the  rela¬ 

tively  small  variation  in  magnitude  of  S^/Sw  compared  with  that  of 
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The  variation  of  CfRe^/Nu  is  from  zero  to  7.4  for  a  surface  of 
temperature  twice  the  free -stream  stagnation  value  and  from  zero  to 
2.8  for  a  surface  held  at  a  temperature  of  absolute  zero,  as  shown  in 
figure  10.  This  indicates  the  inadequacy  of  utilizing  the  flat-plate 
value  of  2.0,  as  has  often  been  done  for  estimates  of  heat  transfer. 
Figure  10  is  of  particular  use  in  evaluating  the  heat  transfer  for  a 
problem  when  used  in  conjunction  with  simple  methods  for  determining 
Cf,  as  proposed,  for  example,  in  reference  1. 


SUMMARY  OF  RESULTS 

From  an  analysis  of  the  laminar  compressible  boundary  layer  based 
on  Stewartson’s  transformation  and  including  effects  of  heat  transfer 
and  pressure  gradient,  the  following  results  were  obtained: 

1.  If  the  condition  of  similarity  is  required  and  the  Prandtl 
number  is  constant  but  different  from  1.0,  the  external  Mach  number 
must  be  either  zero,  constant,  or  very  large.  If  the  Prandtl  number 
is  taken  as  1.0,  the  Mach  number  may  be  arbitrary.  The  free-stream 
velocity  distributions  consistent  with  the  similarity  concept  are 
either  power-law  or  exponential  distributions  in  the  transformed  coor¬ 
dinates.  Since  the  exponential  distribution  appears  to  be  limited  to 
favorable  gradients  and  in  this  range  the  problem  may  be  reduced  to  a 
special  case  of  the  power-law  distribution,  the  calculations  have  been 
based  on  the  latter  class. 

2.  For  flows  with  favorable  pressure  gradients,  unique  solutions 
were  obtained.  For  flows  with  adverse  pressure  gradients,  two  types 
of  solution  were  obtained  which  have  been  identified  as  either  essen¬ 
tially  viscid  or  inviscid  in  the  outer  portions  of  the  boundary  layer. 
The  inviscid  solution  sometimes  involved  velocity  overshoot  within  the 
boundary  layer.  For  favorable  pressure  gradients,  the  viscid  solution 
is  required  by  the  boundary  conditions.  For  adverse  pressure  gradients 
there  are  two  viscid  solutions;  these  correspond  to  the  maximum  and 
minimum  wall  shear,  which  exclude  velocity  overshoot. 

3.  For  heated  surfaces  with  favorable  pressure  gradients  a  veloc¬ 
ity  overshoot,  which  increases  with  increasingly  favorable  gradient,  ' 
results  within  the  boundary  layer.  This  excess  velocity  is  associated 
with  the  acceleration  of  a  layer  of  fluid  in  the  outer  portion  of  the 
boundary  layer,  with  density  less  than  the  external  density.  Since 
this  layer  is  subject  to  the  external  pressure  field  and  is  restrained 
only  slightly  by  the  viscous  forces  acting  on  it,  it  is  accelerated 
more  than  the  external  flow. 

4.  For  a  Prandtl  number  of  1.0,  when  the  thicknesses  of  the  dy¬ 
namic  and  thermal  boundary  layers  are  defined  by  a  fixed  fraction 
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(say  0.99)  of  the  velocity  ratio  or  stagnation-temperature-difference 
ratio,  the  thermal  boundary  layer  is  thicker  than  the  dynamic  layer 
for  favorable  pressure  gradients  and  thinner  for  adverse  gradients. 

5.  The  variation  of  a  Reynolds’  analogy  parameter  is  from  zero  to 
7.4  for  a  surface  of  temperature  twice  the  free-stream  stagnation  value 
and  from  zero  to  2.8  for  a  surface  held  at  a  temperature  of  absolute 
zero,  with  the  value  2.0  for  the  flat  plate. 


Lewis  Flight  Propulsion  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Cleveland,  Ohio,  October  15,  1954 
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APPENDIX  A 


SYMBOLS 


The  following  symbols  are  used  in  this  report 


sonic  velocity 


C,Ci,C2,  etc.  arbitrary  constants 


local  skin-friction  coefficient ,  Cf  = 


specific  heat  at  constant  pressure 
function  related  to  stream  function  by  f  = 


asymptotic  function,  g  =  f^ 
enthalpy 


thermal  conductivity 


Sutherland  r  s  constant 


local  external  Mach  number,  M( 
exponent  from  Ue  =  CX111 


Nusselt  number,  Nu  = 


Prandtl  number,  Pr  = 


*0  -  Sr 


[  m  +  1 
2V0UeX 


static  pressure 
Reynolds  number,  Re 


pwuex 


enthalpy  function,  S 
static  temperature 


transformed  longitudinal  velocity  component. 
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2  oscillation  coefficient,  eg,  (C2) 

co  damping  coefficient,  eg.  (C3) 

Subscripts : 

e  local  flow  outside  boundary  layer  (external) 

J  result  of  j’k*1  iteration 

s  stagnation  value 

w  wall  or  surface  value 

0  free-stream  stagnation  value 

Other  notations: 

~  asymptotic  guantity 

primes  denote  differentiation  with  respect  to  T] 
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APPENDIX  B 


ASYMPTOTIC  SOLUTION 


To  evaluate  the  integrals  in  equations  (24)  and  (27),  it  is  nec¬ 
essary  to  have  closed-form  expressions  for  the  integrands  concerned,  in 
the  range  of  large  tj.  This  requires  a  Solution  of  the  system 

f"’  +  ff"  =  0(f'2  -  1  -  S)  (18a) 

S"  +  fS'  =0  (18b) 

for  large  tj,  which  is  €he  asymptotic  solution. 


The  asymptotic  solution  for  f  (designated  f)  Is  assumed  to  con¬ 
sist  of  a  sum  of  terms,  each  smaller  than  the  preceding.  Only  the 
first  two  terms  will  be ^discussed  herein.  The  corresponding  solution 
for  the  enthalpy  term  S  is  also  obtained. 


Let 


/V  IV 

f  =  fl  +  f2 


(Bl) 


where 


•V  M 

f2  «  fX 

?2  «  *i 


Now,  since  lim(f')  =  1,  let 

T)-*00 


-  X 


(B2) 


where  x  is  an  undetermined  constant.  If  f-^  is  inserted  into  (18), 
the  corresponding  enthalpy  term  S-,  must  be  identically  zero.  In¬ 
serting  equations  (Bl)  and  (B2)  into  equations  (18)  and  dropping 
higher-order  terms  result  in 


+  (n  -  x)f->  =  p 
S’’  +  (7]  -  X)S*2 


(B3) 


The  energy  equation  can  be  Integrated  directly  to  give 

-itsi! 

Si  -  0»  2 
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which  integrates  once  again  to  the  complementary  error  function 
noted  cerf) 


(de- 


(B4) 


If  equation  (B4)  ia  now  substituted  into  the  momentum  equation  of  equa¬ 
tions  (B3),  with  the  notation 

g(n)  = 

there  results 


g"  +  (t)  -  x)g' 


A  particular  integral  to  equation  (B5)  is 


(B5) 


g  = 


(B6) 


The  complementary  function  can  be  found  by  noting  that  the  homogeneous 
part  of  equation  (B5)  is  Weber's  equation.  Hartree  (ref.  7)  gives  the 
general  solution  for  large  values  of  the  argument  (q  -  x)  which  can  be 
written 


« 


,  . -(2p+l) 

g  =  ^(ti  -  x)  exp 


■to.  Jfr 


+  OgCq  -  x) 


2P 


where  and  c are  undetermined  constants. 


(B7 ) 


For  p  0  it  is  clearly  necessary  to  take  =  0  if  the 
boundary  condition  lim  g  =  0  is  to  he  applied.  For  p  <  0  the 

T)-H» 

boundary  condition  does  not  require  =  0;  this  introduces  a  lack  of 
uniqueness  in  this  range.  The  significance  of  =  0  was  more  fully 
discussed  in  the  section  1OTIQUEM1SS. 


* 
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Using  the  first  term  of  the  expansion  for  the  complementary  error 
function 


and  combining  the  preceding  equations  result  in  the  following 
expressions : 
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APPENDIX  C 


CONVERGENCE  AND  EXTRAPOLATION 


the  os- 


place  of  f!|+1  to  obtain  f^, 

dilation  ia  reduced  and  a  convergence 
takes  place.  A  typical  result  is  shown 
in  sketch  (b). 

When  the  value  for  (3  for  which  a 
solution  was  sought  was  sufficiently  posi¬ 
tive  ,  the  enthalpy  function  S  also  showed 
a  tendency  to  oscillate.  In  these  cases, 
applying  the  same  averaging  procedure  to  S 
again  improved  the  convergence.  It  was  also 
found  that  convergence  was  improved  if ,  in 
the  intial  assumed  function  for  f"(f'),  the 


Iteration 


The  method  of  successive  approximations  used  in  solving  equations 
(24)  and  (28)  is  as  follows:  Two  functions  fj(f')  and  Sj(fr)  are 

assumed  and  inserted  into  the  right  sides  of  equations  (24)  and  (28). 
This  produces  two  new  functions,  fj+]_(f ')  a^id  Sj+-^(f ')  on  the  left. 

The  question  of  convergence  is  the  first 
to  consider.  In  reference  2,  Crocco 
treated  a  momentum  equation  which  was 
essentially  equation  (24)  with  (3  =  0. 

There  it  was  shown  that  the  result  might 
converge  to  a  pair  of  functions  between 
which  it  would  oscillate  and  of  which 
the  geometric  mean  was  the  proper  solu¬ 
tion.  In  practice,  the  use  of  the  arith¬ 
metic  mean  was  demonstrated  to  be  adequate. 

In  the  same  way  in  the  present  case,  the 
property  of  oscillation  cannot  be  developed 
analytically;  however,  it  has  been  found 

f"  +  f" 

by  trial  that,  if  J  J+A.  is  used  in 


f ' 
(b) 


3lope  was  taken  so  that  it  satisfied  equation  (18a);  that  is, 


©. 


-P(l  +  sv) 

f  " 
w 


When  an  iterative  method  is  used  to  determine  a  function,  it  is 
always  desirable  to  develop  a  method  of  extrapolating  the  result  to 
correspond  to  a  larger  number  of  iterations  than  have  actually  been 
carried  out.  This  cannot  be  done  in  an  exact  fashion  unless  a  definite 
law  of  convergence  Is  established.  Recently,  an  extrapolation  method  was 
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devised  (ref.  24)  which  required  four  successive  iterants  for  an  arbi¬ 
trary  iterative  computing  scheme.  The  development  assumed  that  the 
remaining  error  after  any  iteration  consisted  essentially  of  two  terms , 
both  of  which  damped  by  a  factor  co  with  each  iteration.  The  sign  of 
one  of  these  terns  was  assumed  to  change  with  each  iteration.  This 
method  extrapolated  a  function  by  breaking  it  into  n-1  parts  and 
treating  it  somewhat  like  an  n-dimensional  vector.  The  method  has 
been  demonstrated  for  Laplace rs  equation  for  which  it  was  quite  adequate. 
For  nonlinear  equations,  however,  the  method  is  not  as  suitable. 

In  referep.ce  1,  a  method  requiring  five  successive  iterants  was 
developed  which  combined  the  method  of  reference  24  and  the  geometric 
mean  rule.  The  function  to  be  extrapolated  is  considered  to  be  made 
up  of  a  set  of  numbers  F^,  where  the  subscript  i  identifies  the 
particular  component  of  the  set.  Then,  the  resulting  relations  for  the 
i^h  component  of  the  extrapolated  function  F  in  terms  of  the  pre¬ 
ceding  five  iterants,  (F*)  .  •••  (Fi)j+4,  where  j  is  the  iteration 
number,  are :  J 


where  the  oscillation  coefficient  9 .  is  given  by 


and  the  damping  coefficient  oo  is 
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Application  of  this  system  was  extremely  effective.  It  generally 
reduced  the  oscillation  remaining  after  five  iterations  by  a  factor  of 
10.  A  typical  plot  of  the  oscillation  of  f^  is  indicated  in 
sketch  (c). 
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APPENDIX  D 


CALCULATION  PROCEDURE 

The  successive  approximation  calculations  were  carried  out  by 
means  of  IBM  Type  604  Calculating  Punch  machines.  The  program  was 
coded  for  fixed-point  calculation,  with  the  standard  Function- 
Generating  control  panel  used,  plus  a  control  panel  especially  wired 
for  rapid  integration  of  quotients  hy  a  trapezoidal  rule.  The  step  size 
(in  f ')  varied  from  &  maximum  value  of  0.050  or  0.025  at  fr  =  0  to 
0.00001  at  f T  =  0.9999,  the  total  number  of  intervals  being  122  in 
the  former  case  and  236  in  the  latter.  By  doubling  and  halving  the 
step  size  for  a  critical  case,  the  results  are  judged  to  contain  a 
maximum  error  of  0.0002.  Comparison  with  solutions  obtained  by  for¬ 
ward  integration,  for  the  same  case,  confirms  this  accuracy.  A  given 
iteration  (utilizing  the  0.050  step  size)  could  be  carried  out  in 

approximately  1—  hours  by  an  experienced  machine  operator.  If  the 

averaging  and  extrapolation  techniques  described  in  appendix  C  are  used, 
10  iterations  generally  would  suffice  for  the  accuracy  desired.  In 
contrast  with  forward  integration,  this  number  of  iterations  is  not  a 
function  of  the  experience  of  the  person  carrying  out  the  calculations. 

In  the  derivation  of  the  integral  relations  (eqs.  (24)  and  (27)), 
it  was  assumed  that  the  velocity  ratio  varied  smoothly  and  monotoni- 
cally  from  zero  at  the  wall  to  1.0  at  infinity.  However,  in  the  range 
P  >  0  and  >  0  (favorable-  pressure  gradient  and  hot  wall),  the 
solution  involves  an  increasing  velocity  ratio  to  a  value  greater  than 
1.0,  followed  by  a  smooth  decrease  to  1.0.  Under  these  unusual  cir¬ 
cumstances,  the  method  of  successive  approximation  derived  herein  must 
be  considerably  modified  if  it  is  to  be  used  at  all.  For  these  cases, 
forward  integrations  were  performed  by  Dr.  Lynn  U.  Albers. 

Equations  (18),  ‘together  with  the  boundary  conditions  (17),  consti¬ 
tute  a  nonlinear  two-point  boundary- value  problem.  Cases  of  this 
boundary- value  problem  were  solved  by  forward  integration,  with  the 
IBM  Card -Programmed  Electronic  Calculator  (CPC)  used  to  integrate  with 
five-point  integration  formulas. 

For  the  cases  where  the  solutions  are  not  unique  (p  <  0),  the 
solutions  were  obtained  in  two  patterns:  In  one  pattern,  p  and  Sw 
were  fixed  and,  for  a  set  of  values  of  f^,  the  quantity  was 

altered  until  boundary  conditions  at  infinity  were  apparently  satis¬ 
fied.  In  the  other  pattern,  f^  and  S^.  were  fixed  and,  for  a  set 
of  values  of  negative  p,  the  quantity  was  altered  until  boundary 

conditions  at  infinity  were  apparently  satisfied.  An  attempt  was  made 
with  both  patterns  to  include  the  solution  with  the  minimum  value 
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of  the  maximum  velocity  ratio  f^ax  within  the  boundary  layer.  Ex¬ 
cept  for  those  cases  where  no  solution  existed  without  velocity  over¬ 
shoot,  this  minimum  value  was  1.0. 

The  details  of  the  integration  method  used  are  described  very 
completely  by  Lynn  U.  Albers  in  an  appendix  to  reference  25.  The 
possible  error  contained  in  the  results  is  indicated  in  the  footnote 
to  table  I.  Each  trial  run  of  a  case  required  approximately  30  minutes. 
A  person  considerably  experienced  with  the  method  of  obtaining  solu¬ 
tions  by  forward  integration  generally  achieved  convergence  within  12 
trials;  however,  tests  indicate  that  this  number  is  insufficient  by  a 
factor  of  the  order  of  2  if  the  person  lacks  experience. 
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tc  be'  correct  to  +0.0002.  Where  f”  and  Sw  were  initially  fixed,  0  and  S^  are  believed  to  be  correct  to  ±0.0002  (except 
in  the  case  of  0  =  0.2460,  Sw  =  -0.4,  where  0  and  S^  are  believed  to  be  correct  to  ±0.002).  The  values  in  the  tables 
are  of  comparable  accuracy  except  at  large  T},  where  the  entries  may  contain  errors  as  large  as  twice  the  above  amounts. 
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-  .0  1  0  9 

.0  09  4 

-  .0  2  6  3 

2.9 

2.6  5  8  7 

1.0  0  3  2 

-.0  0  8  6 

.0  07  0 

-  .0  2  0  2 

3.0 

2.7  5  8  9 

1.0  0  2  4 

-.0  0  6  6 

.0  05  3 

-.0  15  4 

3.1 

2.8  5  9  1 

1.0  0  1  8 

-.0  0  5  1 

.0  039 

-.0  117 

3.2 

2.9  5  9  3 

1.0  0  1  4 

-  .0  0  3  8 

.0  02  9 

-.0  0  8  7 

3.3 

3.0  5  9  4 

1.0  0  1  1 

-  .0  0  2  8 

,0  022 

-.0  0  6  4 

3.4 

3.1  5  9  5 

1.0  0  0  8 

-.0  0  2  1 

.0016 

-.0  0  4  7 

3o5 

3.2  5  9  6 

1.0  0  0  6 

-  .0  0  1  5 

.0012 

-.0  0  3  4 

3.6 

3.3  5  9  6 

1.0  0  0  5 

-  .0  0  1  1 

.0  00  9 

-.0  0  2  5 

3.7 

3.4  5  9  7 

1.0  00  4 

-  .0  0  0  7 

.0  00  7 

-.0  0  18 

3.8 

3.5  5  9  7 

loO  0  0  4 

-  .0  0  0  5 

.0  00  6 

-.0  0  12 

3.9 

3.6  5  9  8 

1.0  0  0  3 

-.0  0  0  3 

.0  00  4 

".0  0  0  9 

4.0 

3.7  5  9  8 

1.0  0  0  3 

-  .0  0  0  2 

.0  00  4 

-.0  0  0  6 
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TABLE  II.  -  SUMMARY  OF  HEAT -TRANSFER 


AND  WALL-SHEAR  PARAMETERS. 


Sw 

e 

fM 

w 

S  1 
w 

C.Re 
f  w 

Nu 

-1.0 

-0.326 

0 

0.2477 

0 

-.3657 

.050 

.2958 

.3381 

-.3884 

.140 

.3527 

.7939 

-.360 

.2448 

.4001 

1.224 

-.30 

.3182 

.4262 

1.493 

-.14 

.4166 

.4554 

1.830 

0 

.470 

.470 

2.000 

.50 

.5806 

.4948 

2.347 

2.00 

.7381 

.5203 

2.837 

00 

o 

1 

-0.10 

-0.0686 

0.0447 

-2.456 

-.2685 

-.050 

.1829 

-.4374 

-.3088 

0 

.2261 

0 

-.325 

.0493 

.2545 

.3100 

-.3285 

.0693 

.2644 

.4194 

-.3285 

.110 

.2818 

.6245 

-.325 

.1354 

.2913 

.7438 

-.30 

.2086 

.3155 

1.058 

-.14 

.3841 

.359 

1.712 

0 

.470 

.376 

2.000 

.50 

.6547 

.403 

2.599 

1.50 

.8689 

.4261 

3.263 

2.00 

.9480 

.4331 

3.502 

-0.4 

-0.235 

-0.050 

0.0447 

-0.8949 

-.246 

0 

.1249 

0 

-.2483 

.050 

.1360 

.2941 

-.24 

.1064 

.1474 

.5775 

-.20 

.2183 

.1626 

1.074 

0 

.470 

.188 

2.000 

.50 

.7947 

.209 

3.042 

2.00 

1.3329 

.2304 

4.628 

0 

-0.1947 

-0.050 

0 

a 

-.1988 

0 

0 

0 

-.16 

.1905 

0 

.9480 

0 

.470 

0 

2.000 

.50 

.9277 

0 

3.436 

1.00 

1.2326 

0 

4.317 

1.60 

1.5213 

0 

5.122 

2.00 

1.6870 

0 

5.565 

1.0 

-0.10 

-0.1613 

-0.2076 

-1.554 

-.1305 

-.050 

-.3139 

-.3186 

-.1295 

0 

-.3388 

0 

-.10 

.1805 

-.4033 

.8956 

0 

.470 

-.470 

2.000 

.30 

.9829 

-.5457 

3.602 

.50 

1.2351 

-.5725 

4.315 

1.00 

1.7368 

-.6154 

5.644 

1.50 

2.1402 

-.6425 

6.662 

2.00 

2.4878 

-.661 

7.527 

aThis  value  was  not  calculated. 
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Figure  1.  -  Family  of  solutions  for  adverse  pressure  gradient.  Sw  =  -0.8. 


sinim 

Kssssssssssttsssssss 


■BBBBBKBBI 

ihbbbhbcibbi 


S^SSSSSHSSSSBSSSIBS 

IIB^BB^BKBBBB^BBBBBBBBBH 

YBBft£BB£l*3aBfl£BBEiBBBBBfl] 

UVBBtfilSBBKBB^BBB^BBBBBBBa 

lYIMBMtaBBBBK'lBBttBBaBBBil 

l^^flBB^B^BB^BB^BB^BBBBBil 

B?.B3IBaB£B£BBBBB^BB^BBBBI 

B»^ISBBBKSIK9IBH!UMKBBBBI 

BBBBBBaBBBBaaBBaBBBftBBBI 


VBBBj 

fi.iiBBgj 
iVBBj 
IHB 

EktUflfl 

■  tHVIBI 

ibbbb^bsbM 
UHiMKvinnviinn 

.....  ■UVBKMmBBBRflfl 
liiBBBBiRBKBSlBBaBBBBBi 

BBBB  BBBBBB9ISBBBBBBBC 
BBBI  BBBBSRBBiSBiaiBBBBBiiHH 
BBBiIBBaBBgBBiiaaBBiBBBBBa 

MmmtMumnumxmmmmmmamamam 

Bbbuiibbbsibbbbbbbbjbbbbbb 

IIMMflBMBMMBHH 

liBBBBJKBBBIgBBBBBBBigliBBBBl 

BBBBBB3B(IBBBSiBBB3BaB*aM 
BBB.lBRB^BBBBBfeSBBflfciBBKBBB] 

BBBKBBBBKSBBBBBfiBBBSBKVBl 

BBBBBSBBBBIiSaBBBBSaflBBk'Uj 

|BBBBBBBggBBBBggBBBi>aBlgBJM| 

■■■antti!!  BB  £5  SB  K9  IS  VI 

l5BBBBBBBBBgBa£g*g^l 


IV1BBBBI 

■VBBBBI 


IBBBBBBBg£S 


*  *  H 


\wmm 

m 


HBK1M 

UgBKBSi 

IRbbbbbbi 

UiraBBBBBISBBLBBBBBBBBBI 

^BMfeSBBBnB^B^BBBBBBBI 

fissssssssss^^ssssi 

lIBBBBgEBlgaBglBgiSlM 

KjpgSiSHB 

kgBBKgBBBgBBiSBBBklill 


sssssss;ss3s:ss:ssi;s 

BBBIKSBBaBKaiSBBfilBKBn 


CBV&  back 


III 

IIBBI 
11BBI 

VBBI 

iflij 

mi 

nil 

iWkwB 

mwSm 

ivuaBiSBKSSBwaSBiiiiw 
[xvkxwzKmaxmwmmmmmwm 
ImvAxizamKititxmmmmnmM 
mv-»Z4n*’*Km  aBBBBBBBBB 

Inanni  ibbbi 

■NCBVtfKBBI 
[BKHKttKKBZil 


IBBSBSKBMI  IBBBI 

WiiiBBK^KSSSBBI  UM 

- 

IBBB^^sXXll 

bbbbb^^Ib! 

^■bbhhb"! 


—■M 

IbbbbbbbbbbbbbbbbbbbI 


III 

mill 

IVBI 

llfll 

tin 

IVKBI 

nanfl 
\mmm\ 

MBBKI 
HI 

■Kami 
p\i«ni 


I^Mbbbbi 

mSKiiii 

i^sanam  n 

IBB^SSKSim  IBB 

— BBBBB^SSSagBBBB— 
IBBBBBBgSSglBMBBI 

■BngnHHn&^s^BBS 

IBBi^Snj 

|BBBBh.^j 


60 


MCA  TN  3325 


g  'uoTqounj  XrIxT3uq.ua 


Figure  4.  -  Enthalpy  function  representation  in  velocity  plane.  (*  denotes  lower- branch  solutions. 
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Pressure- gradient  parameter,  (3 


(a)  Favorable  and  adverse  pressure  gradients. 
Figure  6.  -  Effect  of  pressure  gradient  on  rail  shear 
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Figure  7.  -  Stagnation  enthalpy  gradient  across  boundary  layer.  (*  denotes  lower- branch  solutions. 


Stagnation  enthalpy  gradient  at  wall 
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Figure  8.  -  Variation  of  heat  transfer  with  pressure  gradient. 
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